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Abstract 

We propose extension of the numerical method to model effect of Bose- 
Einstein correlations (BEC) observed in hadronization processes which allows 
for calculations not only correlation functions C2{Qinv) (one-dimensional) but 
also corresponding to them C2(Q x ,y,z) (he., three-dimensional). The method 
is based on the bunching of identical bosonic particles in elementary emitting 
cells (EEC) in phase space in manner leading to proper Bose-Einstein form 
of distribution of energy (this was enough to calculate C^Qinu))- To obtain 
also C2(Q x ,y,z) one has to add to it also symmetrization of the multiparticle 
wave function to properly correlate space-time locations of produced particles 
with their energy-momentum characteristics. 

Key words: hadronization • correlations • BEC 

Recently we have proposed novel numerical method to model effect of the Bose- 
Einstein correlations (BEC) observed in hadronization processes (see |2] for de- 
tails and for other references concerning generalities of BEC). This method allows 
for calculation of the so called invariant (i.e., one-dimensional) correlation function 
C2{Qinv)- Referring to P3I2] for motivation and justification and all other details let 
us list here only the main points proposed so far: 

• Bosonic character of produced secondaries demands that they are produced 
in bunches (named by us Elementary Emitting Cells - EEC's) in momentum 
space jS]. 

• The problem then is how to model production of EEC. It is done by adding 
to some preselected particle with energy E\ (according to some distribution 
f(E), which in our numerical calculations was taken in the Bolztmann form, 
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f(E) ~ exp(— E/T), with temperature T being a parameter) another particles 
of the same energy with some probability P = Pq ■ exp(—Ei/T) (where Pq is 
another parameter) up to first failure. After it one starts to build another 
EEC and continues as long as total energy allows. 

• Such procedure results in a number of EEC's (distributed according to Poisson 
distribution and follow boltzmanian energy level distribution) with a number of 
identical particles in each of them distributed according to geometrical (Bose- 
Einstein) distribution with energy level distribution following Bose-Einstein 
statistics. 

• In this way one gets C2(Qinv) = 2 but only in the first bin, i.e., for Qi nv = 
0. In order to get the characteristic shape of C 2 (Qi n v) function one has to 
allow for particles in each EEC to have slightly different energies, for example 
distributed around E± according with gaussian distribution with the width a, 
which is our next and last parameter on this stage. It should be stressed at this 
point that in the field theory approach to BEC, as for example that presented 
in a — corresponds to infinite source and commutation relations with 
Dirac-delta functions, whereas nonzero values of sigma arise for finite space- 
time extensions of the hadronizing sources. 

To obtain also correlation functions C2(Q x ,y,z) ( i-e., three-dimensional, calcu- 
lated for different components of the differences of particle momenta Q)one has to 
proceed further. What we propose here is the following: 

• Change of Qi nv to Q means that we shall be now sensitive not only to overall 

spatial difference between particle production points r but to the whole vector 

r as well. We have to then assume that particles are produced from some 

spatial region and that density of production points is given by some function 

(our additional input) p(x, y, z) 1 
^or a time being we are assuming for simplicity that hadronization is instantaneous so p is 
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• Momenta of each particle in a given EEC, p i} obtained in the first part of 
algorithm must now be decomposed into their components, {p^\ Py\ P^) ■ To 
do this let us observe that BE statistics demands that multiparticle wave func- 
tions must be symmetrized accordingly and this results in correlations between 
production points represented by p(x,y,z) and momenta {Px\Py\p z ^)- -De- 
noting by 5i =X! y tZ the corresponding differences in position this correlation is 
given (in the plane wave approximation) by the known 1 + cos (5i ■ Qi) term. 

• We proceed then in the following way. In each EEC we select the (pW , pW ; pM^j 
for the first, i — 1 particle in some prescribe way (here isotropically but one 
can introduce at this moment px cutting or something else as well) and then 
establish (p^\p y l \p z 1 ^ f° r every one of additional particles, i > 2, in such a 
way that cos(<5j • Qi) < 2 • Rand — 1 where Rand is random number uniformly 
chosen from interval (0, 1). 

This leads us to results presented in Figs. 1 and 2. We regard them as very promising 
but we are aware of the fact that our proposition is still far from being complete. 
To start with one should allow for time depending emission by including 5E ■ St 
term in the cos(. . .) above. The other is the problem of Coulomb and other final 
state interactions. Their inclusion is possible by using some distorted wave function 
instead of the plane waves used here. Finally, so far only two particle symmetrization 
effects have been accounted for: in a given EEC all particles are symmetrized with 
the particle number 1 being its seed, they are not symmetrized between themselves. 
To account for this one would have to add other terms in addition to the cos(. . .) 
used above - this, however, would result in dramatic increase of the calculational 
time. 

time-independent . 



4 



c 



©; 



2.0 

1-9 

1.8 

1-7- 

1.6 

1-5 

1.4 

1.3 

1.2 
1.1 
l.O 

0.9 



T=3.5 GeV, <r = o.3*T GeV, P=o.7*exp(...) 

■ ■ ■ - W= 46.5 GeV 
- w= 91.2 GeV 
W - 182.4 GeV 



0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 



2.0 
1.9 
1.8 

1-7 
1.6 
1.5 
1-4 
1-3 
1.2 
1.1 
1.0 



0.9 

0.0 0.2 



•T = 3.1 GeV 
T - 3.5 GeV > a - o,3*T GeV, P=o.7*exp(..0 
T- 3 . 9 GcV 




0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 
Q. [GeV] 



- - - - p o = 0.6 1 

■ P° = 0.7 > T = 3.5 GeV, a=o,3*T GeV 

P° = 0.8 I 




0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 

Qi„v ^eV] 



2.0 
1-9 
1.8 
1-7 
1.6 
1-5 + 
1.4 
1.3 
1.2 
1.1 
l.O 
0.9 



" CT, = O.l 

1° = 0.3 > T = 3.5 GeV, P=o.7*exp(..J 
"„ = 0.5 I 




0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 

Q. [GeV] 



S 11 




-2.0 -1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 2.0 



Q [GeV] 



el 




2.0 -1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 2.0 
Q [GeV] 



o; 




-2.0 -1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 2.0 

Q xv2 [GeV] 



C 1.1 



;.o -1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 2.0 

<5 [GeV] 



Fig.l Example of results obtained for C2{Qinv) (left column) and corresponding C2(Qx,y,z) (right column, for 
parameters used here all C'2{Qx,y,z) are the same). Calculations were performed assuming spherical source p(r) of 
radius R = 1 fm (spontaneous decay was assumed, therefore there is no time dependence) and spherically 

symmetric distribution of p x ,y,z components of momenta of secondaries p. Energies were selected from 
f(E) ~ exp(— E/T) distribution. The changes investigated are - from top to bottom: different energies of 
hadronizing sources, different temperatures T, different values of parameter Po and different spreads a = ctq ■ T of 

the energy in EEC. 
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Fig. 2 Demonstration of great role played by parameter a denning the energy spread of particles in the EEC. 



To summarize - we propose new method of numerical modelling of hadronization 
events in such way as to respect the bosonic character of produced secondaries and 
therefore leading to BEC. It seems to converge in some sense to the proposition 
presented long time ago in [3] , which was, however, in practice impossible to be im- 
plemented. The only hope that it could work now is that in our case symmetrization 
is within given EEC and not for the whole bunch of particles produced. Therefore 
the number of terms involved is rather limited, whereas in [3] the whole source 
had to be symmetrized at once. But the effect of including at least terms when 
symmetrization between, say, particles 2 and 3 are added to the already present 
symmetrization between 1 and 2 and 1 and 3, must be carefully investigated and 
estimated before any conclusion is to be reached. 
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